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We extent the previous determinations of nonasymptotic critical behavior of Phys. Rev B32, 7209 
(1985) and B35, 3585 (1987) to accurate expressions of the complete classical-to-critical crossover 
(in the 3-d field theory) in terms of the temperature-like scaling held (i.e., along the critical isochore) 
for : 1) the correlation length, the susceptibility and the specihc heat in the homogeneous phase for 
the n-vector model (n = 1 to 3) and 2) for the spontaneous magnetization (coexistence curve), the 
susceptibility and the specihc heat in the inhomogeneous phase for the Ising model (n = 1). The 
present calculations include the seventh loop order of Murray and Nickel (1991) and closely account 
for the up-to-date estimates of universal asymptotic critical quantities (exponents and amplitude 
combinations) provided by Guida and Zinn-Justin [J. Phys. A31, 8103 (1998)]. 

05.70.Jk, ll.lO.Hi 


I. INTRODUCTION 


The asymptotic critical behavior, characterized by universal quantities (exponents and amplitude combinations), 
is now theoretically well established with accuracy However, the comparison of the theoretical results with 
experimental or numerical data is made easier when the theoretical expressions are extended into regimes where 
the asymptotic pure scaling breaks down iH (calculations done far away from the critical point, characterized by 
nonuniversalities and including eventually crossover phenomena, see reviews p|J^-|ll||). This extension has appeared 
necessary notably when measurements on colloids [Q, but also on complex systems such as ionic fluids [Q or 
polymers [Q, seemed to yield strong nonuniversalities in approaching the critical point. Indeed theoretical studies 
have suggested that, in some cases, those nonuniversalities could be due to a phenomenon of “retarded criticality” 
which characterizes measurements done outside the asymptotic critical domain |p^ , |l^ . Several recent theoretical 
(and/or numerical) studies have as well explicitly considered the evolution of effective exponents with emphasis 


on their monotonic or non-monotonic character 
crossover for Ising systems is not yet clear-cut 
nonasymptotic critical behavior from held theory 


11|J17|-|21|. Furthermore the description of the classical-to-critical 
. For these reasons and because our preceding determination of 
did not yield continuous functions covering an entire crossover 


region, it seems to us useful to consider again those calculations in order to (see also ^E&y- 


1. extend them to a complete account of the classical-to-critical crossover which characterizes the framework of 
Held theory [ p6p^ ], 

2. include the seventh order series for the critical exponents determined by Murray and Nickel in order to 
account as closely as possible for the up-to-date estimates of universal asymptotic critical quantities (exponents 
and amplitude combinations) provided by Guida and Zinn-Justin (referred to in the following as GZ). 

In the previous work and contrary to an initial attempt regarding the homogeneous phase (n = 1), we 

provided only continuous expressions of t valid for t ^ 10“^ {t is the temperature-like scaling Held which is proportional 
to the absolute value of the reduced critical temperature |T — Tc\ /Tc). The crossover was not completely described 
because it was thought at that time that the Held theoretical framework had a range of validity strictly limited to the 
Hrst correction to scaling term. Gonsequently, the practical limit of physical validity of the functions was imposed 
by the range of t where the second correction to scaling term speciHc to Held theory becomes non negligible and this 
occurs about t ~ 10“^. Since then, it has appeared that the range of validity of Held theory could be much larger 
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than that and even could cover the entire classical-to-critical crossover region ||T^ that it describes. Somehow it is 
interesting to give expressions valid in the entire crossover region if only because it may be compared to other kinds 
of classical-to-critical crossovers either experimental |p^ or from numerical studies which, under some 

particular conditions, are identical |ll8| , |l9 , ^4|J^J^ to the field theoretical form (but see also our comment in 

For technical reasons, in the preceding work we did not constrain our theoretical expressions to include 
very closely the estimates of universal asymptotic critical quantities of that time (with their error bars) so that 
uncertainties were underestimated and then our estimates of the correction amplitude ratios were, presumable also 
not firmly determined. Moreover small errors existed in the preceding study of the inhomogeneous phase (as 
indicated elsewhere [18p2P^) which have been eliminated from the present work (nevertheless, we have explicitly 
verified (see fig. 1 of ref. ||33(]) that the errors have had no important consequence on the final results as it could be 
clearly deduced from a comparison of our estimates of universal amplitude-combinations with those of Guida and 
Zinn-Justin Q,^). 


II. PRINCIPLE OF THE CALCULATIONS 


A. Brief reminder 


As in the preceding work [|^[^ , and using the same resummation method, we have considered the correlation length 
^ (t) (in the homogeneous phase T > Tc for the n-vector model with n =1 to 3), the susceptibility y (f) and the specific 
heat C (t) (in the homogeneous phase T > Tc with n =1 to 3 and in the inhomogeneous phase T < Tc with n =1) 
and the coexistence curve (spontaneous magnetization) M (t) (in the inhomogeneous phase with n =1). 

For practical reasons, it is useful to fix our notations relative to the actual asymptotic critical behaviors (i.e., in 
terms of the physical variable r = ^ 0^ instead of t ^ 0, see also section III): 


CM 

xM 

C(r) 
M ij) 




l + a±|r|^+0(|rr) 


= C 

= r± |r|-^ [l + a± |r|^+0(|rp^)] 

= — [1 + aa^ \t\^ + O (|Tp^)] -b Be 


= B\t[ 


1 + om ItI -b O (|rM) 


(la) 

(lb) 

(l c) 

(l d) 


in which a, /3, 7 and iz are the critical exponents, A (also denoted hy 9 = ujv by GZ) is the correction exponent, 
F±, and B are the leading critical amplitudes and o^, a^, and qm are the (confluent) first-correction 
amplitudes, finally Bcr is a critical background. One usually restricts the consideration of the critical singularities 
to small values of t oc |t| as it is implicitly assumed in Eqs. #■ The obtention of nonasymptotic critical behavior 
supposes the explicit consideration of non-necessarily small values of t 

Let us suppose that we want to calculate the susceptibility y as function of the (non-necessarily small) temperature¬ 
like scaling field t. Galculations of such a nonasymptotic critical behavior from (the massive) field theory (in three 
dimensions) 1-0 present the following features (additional details may also be found elsewhere p5|): 


1. the function y (t) is primarily performed under the implicit form because the quantities y and t are primariN 

given as perturbation series in powers of the renormalized coupling parameter g (up to fifth or sixth Q| 

order): the functions y {g) and t (g) are resummed for g varying in the range ]0, g*[ where g* is the zero of the 
Wilson function (or the “/3-function”) W(g) also primarily given as power series of g (up to sixth order ; g* 
is the fixed point value of g). 

2. The consideration of discontinuous values of g is a compelling need of the numerical resummation procedure. 
Consequently fitting an ad hoc function of t to the calculated points eliminates the auxiliary variable g and 
provides us with the final expression of y (t) as the explicit continuous function of t we are looking for in the 
range \t\ G ]0, -boo[. 

3. The actual calculation of the quantities of interest [like y (g)] for values of g close to g*, at which point they 
are singular (due to the critical singularity we expect to closely reproduce), requires expressing them under an 
integral representation like: 


X ( 5 ) = X iVo) exp 



V (x) W (x) 


( 2 ) 
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in which 7 (g) and v {g) are not singular at g* and are primarily given as power series of g (up to seventh order 
@). Especially, 7 {g*) and v (g*) provide the field theoretical estimates of the critical exponents 7 and i'. Only 
the elementary series 7 ( 3 ), v (g) and W (g) are resummed using the sophisticated method mentioned in the 
following step [the value yo is chosen small enough to allow a direct simple summation of the series x ( 2 / 0 )] • 

4. To sum perturbation series like 7 (g), v {g) or W (g) for a given value oi g a, Borel-Leroy transformation is 
used, combined with a conformal mapping. An estimation of the error is deduced from the observation of the 
convergence properties of the series when varying the free parameters of the transformation. This leads us to 
fixe those parameters (resummation criteria) in such a way as to obtain a combination of the error bounds on 
e.g., 7 (g), V (g) and W (g) which gives a kind of envelope for y (g) via two functions Xmax (ff) and Xmin (ff) (and 
similarly for t{g)). 

5. Since the critical singularities are similar in the two phases of the transition, the calculations in the inho¬ 
mogeneous phase (T < Tc) do not require the consideration of new series for the exponents compared to the 
homogeneous phase (T > Tc). Hence the same three series 7 ( 5 ), v {g) and W (g) express the critical singularities 
via integrals similar to that given in (^, only new critical amplitude functions of 0 (hence not singular at g*) 
must be calculated [^,[^ and summed using the transformation mentioned in step^ 


B. Improvements to the previous work and presentation of the results 


1. The fitting procedure 


In the present work, the fitting procedure of step of section IIA is performed in the entire range of values 
of g G ] 0 ,( 7 *[. Consequently the entire classical-to-critical crossover specific to the field theoretical framework is 
completely accounted for by our final functions [see Eqs. (^,||) and tables |-IV[. This is illustrated, for the Ising model 
(n = 1 ), by fig. ( 0 ) which displays the evolution of two effective exponents XeS (t) and (t ), which are defined as 
for example [see Eqs. i)]: 


7eff it) = - 


dint 


(3) 


Fig. (j^) shows the effective exponents (calculated from the crossover functions of tables | and i) which interpolate 
between critical and classical (mean field) values following a form of crossover dictated by the framework of field theory. 
Indeed the (massless or critical, i.e., for t = 0) scalar field theory in three dimensions is defined on a special trajectory 
of the renormalization group (a renormalized trajectory [ p^ [RT]) which takes its origin at the Gaussian fixed point 
(characterized by classical values of the “critical” exponents) and joins the Wilson-Fisher fixed point (where the 
critical exponents take on their critical values according to the universality class considered). Of course, the crossover 
so induced is not universal, it is specific of the framework used. In fact, strictly speaking, only the extreme asymptotic 
moving away from the Wilson-Fisher fixed point induced by small non-zero values of t is universal (critical exponents 
and critical amplitude combinations), even the first-correction amplitude (defined in the close vicinity of the fixed 
point when t is not very small) is not universal. For example, in the present work, a specific definite sign of the first 
correction amplitude is imposed due to the RT chosen, however in actual systems that kind of correction may well be 
of the opposite sign and even absent [ p^ . Fortunately, nonuniversal does not mean necessarily absent in actual critical 
behaviors. It may well occur that actual systems (or mode^ display, more or less partially, the kind of crossover 
calculated from field theory [|l^Jl^ , [l^ , p4| , P^J^ . See section IH for a practical use of the cr ossov er functions. 

Let us now give some technical informations on the fitting procedure of step p| of section IIA that must be a pplied 
two times for each quantity considered because of the two error bounds “max” and “min” (see sections IIA and IIB 2 
for the meaning of these bounds). 

In order to analytically reproduce the functions calculated point by point, we use the following general form |^: 


K 

F{t) = Zt‘^Y[(l + Xit°d)y 


+ Xe 


(4) 


in which K is the maximum number of factors (in a preliminary work we had AT = 3, in the present work K can 
be as large as 5), and with: 


£» (t) = A - 1 


Si'i/t 1 

S2'/t 1 


(5) 
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in which A is the correction exponent. We have adjusted each of the parameters Z, e, {AT^, 1^} (i = 1, • • •, K), Si, S 2 , 
Xq and A so as to fit the discretized evolutions of the quantities considered (^, C, Ms) as continuous functions of 
the temperature-like variable t in the range t G . Of course, there are some external constraints on the 

values of these parameters which facilitate their adjustment: 

1 . The exponents e and A must take on values already known from the resummation of the corresponding elemen¬ 
tary series. 

2. The amplitude Z is easily determined with few points corresponding to very small values of t. 

3. To make it easy to get a close reproduction of the crossover towards the classical behavior when t 00 , there 
are constraints: 

• on Si, so that we have (see, for example, Eqs. (A9, AlO) of S. Caracciolo et al. [^): 

1/2 ( 6 ) 

this leads to: 



• on one of the couple {A^, Til’s by imposing that a known classical behavior is reached in the limit t 00 
then it comes: 


e -b 2 'y ' Yi — Cc (8) 

^ i=l 
K 

Z\{{Xi)^^=Z, (9) 

i=l 

with Be and Zf. the classical values of the critical exponents and amplitude respectively. This leads to the 
constraints for one of the {Ai,Ti}’s: 


= 2 (ec - e) - ^ y, 

i/io 

1 


= 


Zr 


n 


-Yi 


i^io 


( 10 ) 

( 11 ) 


with the classical values Bc = 1 , | or 0 and Zc = 2, 1 , -s/O or Be — Xq for respectively the susceptibility, 

the correlation length, the coexistence curve and the specific heat [Xq is the additive critical part of the 
specific heat and its classical value; Sc = 3 in the inhomogeneous phase (T < Tc), while Be = 0 in the 
homogeneous phase (T > T^)]. 


With the above prescriptions, we have been able to reproduce the original calculated points with a maximum 
(local in t) of relative deviation less than 10 “'^ (in the worst case and for a limited number of functions especially 
in the inhomogeneous phase). However, globally (mean value of the local deviations over the entire range of t), the 
adjustment is much better for all the quantities. 

The results of those adjustments to the discrete calculated points are given in tables f-IV. 

We emphasize that the large number of digits displayed in the tables lays no claim to a better accuracy than in the 
work of GZ, it is simply required to obtain a careful fit of the crossover functions to the discontinuous points primarily 
calculated from the available perturbative series. 
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2. The resummation criteria 


In our preceding work the resummation criteria of step of section ilA which gave the bounds “max” 

and “min” were not chosen so as to closely reproduce the uncertainty of the (at that time up-to-date) estimates 
of universal asymptotic critical quantities (exponents and amplitude combinations). They simply proceeded from a 
primary analysis of the convergences of the elementary series [i.e., 7 (g), 1 ^( 5 ), etc... ] resulting from the (unique) 
resummation technique considered. This makes a notable difference because a given function brings several elements 
into play [see, e.g., Eq. (^)] introducing a possible frustration of the individual resummation criteria. Moreover, when 
one determines the error bar for an individual quantity, one often rounds it up because several resummation methods 
may have been considered yielding answers slightly different from each other. Since the various asymptotic critical 
behaviors of the functions of interest (y (t), etc...) result from the combination of a small number of elementary series 
(namely: 7 (g), v (g), W (g) and a few amplitude functions), the individual criteria were combined in our preceding 
work so as to provide an envelope of the error accounting automatically for correlations (frustrations) between the 
error bounds. This has induced some underestimation of the errors when the universal critical exponents or amplitude 
combinations were (re)-considered from the final expressions of the functions compared to their independent estimates. 

The spirit of the present work is different. We have constrained the resummation criteria of the elementary series 
so as to get as closely as possible the GZ estimates for the universal quantities despite the possible frustrations of the 
error bounds mentioned above. Thus we have taken into account the extensions up to seven loops of the series for 
the critical exponents given by Murray and Nickel |28|. For the reasons indicated just above, and also because the 
error estimates of the amplitude combinations of GZ are deduced from the analysis of the parameter dependence in 
the equation of state [Q (they have not been obtained from the direct analysis of specific series for the quantities 
of interest), we have encountered some difficulties in fixing the resummation criteria for some amplitude series (it 
is likely that GZ have overestimated the error for some quantities). In addition, in doing so and concerning the 
amplitudes we have introduced an imbalance between the error estimates of the two phases. Indeed our criteria are 
adjusted so as to get universal ratios (or combinations) of amplitudes which, structurally in the present work, express 
themselves as series strictly defined in the inhomogeneous phase. On the contrary, the resummation criteria in the 
homogeneous phase (for only one amplitude function) have been fixed without constraint. Gonsequently, the resultiM 
error estimates of the correction amplitudes that we presently obtain are larger than in the previous work of refs, [pp] 
and notably in the inhomogeneous phase case; they are presumably overestimated (see below and part IIIB). 

Table ^ shows our estimates of the critical exponents (resulting from our resummation criteria) compared to the 
GZ estimates. One may observe some very small differences due to the fact that, in th e pre sent work, the scaling 
relations are automatically satisfied for each bound “max” or “min” (see step ^of section II A) while only the central 
values of the GZ exponent estimates satisfy the scaling relations (the apparent errors for 7 , v, [3 have been determined 
independently |p|). Table 0 shows how much the respective estimates meet the scaling relations in both cases. Tables 
VI and VII display the values of the universal combinations of leading critical amplitudes as they are accounted for 
by our crossover functions. The degree of agreement with GZ is graphically illustrated by figs (^, ^). 

From tables 110 one may observe that our bounds on the correction exponent A differ from GZ. This is due to 
the correlation of errors mentioned above. Indeed, we have never considered A as an independent constituent of the 
asymptotic critical behavior. Instead it has been (numerically) deduced from the resummation criteria associated 
with the elementary series v (g) and W (g) because of the definition A = ujiy with: 


UJ = 


dW (g) 


dg 

v=v {g*) 


( 12 ) 

(13) 


The resummation criteria for the elementary series W (g) have been chosen so as to yield estimates on the bounds 
on g* very close to those of GZ (see table VIII). 

Similarly, the present determinations of (and the uncertainties on) the first correction-to-scaling terms displayed in 
tables tX and ^ differ from our preceding work essentially be cause of our systematic account for the up-to-date 
estimates of the leading amplitude combinations (see tables ^ and VII). Notice the likely unrealistic smallness of the 
correction ampl itude om in the case “min”. This confirms our probable overestimation of the error on the correction 
terms (see part IIIB). 

It is worth indicating that the values displayed in tables K and are not obtained from the crossover functions of 
tables J-W by simply using the expression [see Eq. (§]: 


K 

aF=J2x,Y, 
2 = 1 


(14) 
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which would be the right expression if the correction exponent in Eq. (Q) was the actual correction exponent A instead 
of the effective exponent D (t) of Eq. (||). To get the values displayed in tables IX and ^ we have made specific fits 

tjll6 functions of EjC^ /1/1IA -i-Viti TTn-i-Vi TO (■f\ - ^ tti T-onrecio TrolntiC' /“if -h ^ 1 O ^ 

preceding work 


) to the theoretical points with D (t) = A in ranges of values of t < 10 ^ (as in the 


III. PRACTICAL USE OF THE CROSSOVER FUNCTIONS 


As already said above, the structural form of the classical-to-critical crossover that we produce here is not universal, 
it is peculiar to the field theoretical framework which corresponds to having performed a limit (the continuum limit) 
in the renormalization group (RG) theory | |4^ ]. The approximation induces the idea that, strictly speaking, the 
“nonasymptotic” calculations would, in fact, be only valid in the close vicinity of Tc- Hence, we do not expect our 
functions to reproduce the experimental data in the entire range t S ]0,+oo[. However, the width £ of the domain of 
agreement between experiments and field theory is not universal: it could actually be reduced (purely and simply) 
to the strict asymptotic critical region (pure scaling laws) or include exclusively the first correction to scaling, but, 
fortunately it may sometimes be much larger and could even cover the entire crossover region! It is our aim to allow 
experimentalists to determine the width £ of the domain of agreement. Notice that, we do not aim at determining (or 
providing) all the ingredients needed to describe the variety of classical-to-critical crossovers that may be produced 
by actual systems, this would be too difficult (due to the infinite variety of nonuniversal contributions). Simply, we 
think our calculation accurate enough to allow the determination of £ for any system allowing, in some sense, the 
determination of subclasses of universality. 

As already explained |5|-|^] , the comparison of the theoretical functions with experimental data involves a very small 
number of adjustable parameters: 

1 . Non-universal global factors, 

2. The proportionality factor 6 between the temperature like scaling field t and r = (neglecting higher 

analytical contributions in r which may sometimes be non-negligible 0 but are out of the scope of our present 
aim): 


t = 9\T\ 


(15) 


3. Additive regular background terms for the specific heat, 

4. Eventually T^. 

For example the comparison with experimental measurements of the susceptibility x may be made as follows: 


Xo Xth {9 \t\) = Xe^p ihl) (16) 

in which Xexp (ItI) represents the experimental data and Xth (t) our function for one Q] of the two bounds “max” and 
“min”. One generally expects theoretically that xo and 9 take on the same values for the two sets of measurements 
above and below provided that the range of values of r is not too large. The fact that xo must be unchanged is 


p+ _l_ 

the consequence of the universality of the ratio with as r ^ 0^: 


Xexp (|t|) Irf^ (l -b rf Irl"^ -!-•••) 

(17) 

in which Tq and Tj*^ are related to our previous definition [Eq. (|lb|)] as follows: 


r± = xo0-^r± 

(18) 


(19) 


As for the stability of 9 this is because the ratio is universal. The values of the universal amplitude combinations 


included in our calculated functions are given in tables VI, VH, [X and ^ 
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A. Redefinition of the role of 9 


If one introduces 9 literally as in Eq. then the fitted leading critical amplitude involves two adjustable pa¬ 
rameters. This is not very suitable. For practical use, we propose to introduce the adjustable parameters as follows 
[compare with eq (Q)] : 


Xexl>(kl) =X0 


^ + +Xe 


t = 0\T\ 


( 20 ) 

( 21 ) 


in which 6 is no longer involved in the pure scaling part of the critical behavior {Z jrj"’^). So introduced, 0 is a nonuni- 
versal parameter which exclusively controls the magnitude of the corrections to scaling. Hence we can progressively 
adjust the theoretical functions to the data starting from the data close to the critical point with 9 = 0 and then 
introducing more and more data with 0^0 (notice that 0 > 0 ) up to the point where consistency is lost. 

The domain £ of r where the experimental data and the field theory agree may involve correction-to-scaling terms 
higher than the first one and this is why it is interesting to have the theoretical expression under the form of a 
complete classical-to-critical crossover (Q . 

Consistency test: If we consider a supplementary set of measurements like the specific heat above and below T^, 
then, by virtue of universality, one must obtain again the same value for 9 with a good fit in a range of values of r 
similar to that considered with y. For the specific heat, it comes: 


|r|) + Ho(r)=Cexp(|r|) (22) 

in which Bq(t) is an additive non critical (i.e., regular or analytic in r) background and Co a nonuniversal multiplicative 
factor which must be the same in the two phases. 

Let us emphasize that the field theoretical form obtained for Cth (t) involves a specific critical additive background 
term which reproduces the famous “classical jump” of the specific heat [see fig (|4|)]. Of course, the magnitude of this 
jump is not universal but in the case where an actual system would reproduce the entire classical-to-critical crossover 
of field theory, then it should also exhibit this jump (up to the global additive background Bo{t) analytic in r). 

If, in addition to C and y, we also possess coexistence curve data, we would have a stronger constraint since then 
no other adjustable parameter would be required to fit those new data. Indeed in the relation: 


Mo Mth (9 |t|) = Mexp (|t|) 


(23) 

Rc and 0 must 


everything is fixed since Mq is related to Co and yo due to the universal amplitude combination |44 
have the same value whatever the quantity considered. 

If we had simultaneously also access to experimental measurements of the correlation length then the constraint 
would be even stronger since again the theory must agree with the data without new adjustable parameter. 

In order to facilitate the use of the crossover functions displayed in table text-files of Fortran code are provided 


B. Account for the error bounds 

We have accounted for the error estimates by providing two sets (“max” and “min”) of functions. In general the 
accuracy of the experimental measurements are much smaller than in the present theoretical calculation so that it is 
not very important to make a difference between the two sets of functions. One or the other choice would provide 
essentially the same quality of the adjustment in the fitting procedure. 

Sometimes accounting for the difference between the bounds “max” and “min” may have some importance so that 
neither one or the other agrees with the measurements but a mixing of the two would. In such a case we propose to 
introduce the mixing via the introduction of a supplementary adjustable parameter E. 

Let us define a new theoretical function as follows 

Fe it) = [Fmax it)f ■ [Emin (t)]^”^ (24) 

which, regarding the definition of the effective exponents corresponds to the linear weighting: 

ceff it)=E- e“// it) + (1 - F) • eT}] it) (25) 
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then the introduction of the other adjustable parameters, such as 9 for example, within Fe (t) is unchanged compare 
to the descriptio n given above. 

As said in part IIB 2 it is likely that the close account of the GZ estimates has led us to overestimate the uncertainty 
on the correction terms so that it seems to us useful to provide also the reader with functions reproducing the complete 
classical-to-critical crossover according to the resummation criteria of the previous (but corrected, see Q) work of 
refs ^,0 although (or rather because), this time, the error are underestimated. This is why we provide two additional 
text-files of Fortran code [ |5| corresponding to the former resummation criteria applied to the corrected series (without 
the seventh order of ref. [ ^ |). 

Acknowledgements We thank Prof. M. Barmatz for interesting correspondence while the present work was in 
progress. 
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FIG. 1. Respective evolutions (calculated from the crossover functions of tables^ and ^ of the effective exponents 7e//(t) 
and aeff{t) in the two phases: the homogeneous (continuous line) and the inhomogeneous (dashed line) phases. Notice, in this 
latter case, the moving down of 7e//(t) below the classical value (= 1.0) in the regime of high values of t. This nonmonotonic 
feature of 7e//(t) in the inhomogeneous phase is in agreement with Refs. |lq,|l9(| and has been numerically observed in Refs. 
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FIG. 2. The two bounds “max” (continuous lines) and “min” (short-dashed line) have been determined so as to reproduce 
as closely as possible the GZ estimates. This is illustrated here with the effective exponents calculated from the crossover 
functions for n = 1 determined in the present work (see tables | and For each exponent, a partial magnification (A) of the 
critical region is provided to show the agreement with the GZ estimates. A similar partial magnification (B) is also provided 
to show the difference with our preceding work of (long-dashed lines). (For other values of n, see table ^.) 
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FIG. 3. The two bounds “max” (continuous line) and “min” (short-dashed line), determined so as to reproduce as closely 
as possible the GZ estimates, account also for the universal amplitude combinations. This is illustrated here with the ratio 


F’^'/r [see Eq. ([lb|)] calculated from the crossover functions for n = 1 determined in the present work (see tables |[l| and 
0). The illustration could have been made as well with the two other universal quantities of table [Vl[ 
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TABLE I. Numerical values of the parameters of the generic crossover function F{t) [Eqs. corresponding to the three 

quantities calculated in the homogeneous phase (T > Tc and for n = 1): the correlation length the susceptibility x S'lid 
the specific heat C. Eor each parameter two values are provided which correspond to the bounds “majc” (upper line) and 
“min” (lower line) of the error treatment. These bounds have been determined so as to reproduce as closely as possible the 
error estimates of GZ on the asymptotic universal quantities (see tables ^ and 0 and text for more details) and provide 
two exclusive sets of functions Emax and Emin- The first row of the table displays the estimate of the universal value of the 
correction exponent A common to all the quantities for the respective bounds “max” and “min”. The values of the parameters 
have been determined by a careful adjustment of F{t) to the discrete evolution of the respective quantities primarily calculated 
by resummation of perturbation series using a Borel-Leroy + conformal mapping method. The specific notations of the two 
leading parameters e (universal critical exponent) and Z (leading critical amplitude) are recalled for each quantity [see Eq. ( 0 )]- 
The symbol “—” means that the term is absent. 



n — 1, homogeneous phase: A 

“max”: 0.49862 
“min”: 0.50516 





C 


e 

P < 

[ 0.631678 

1 0.6290975 

7 \ 

[ 1.2408875 
[1.23830 

—a < 

' -0.1049675 
-0.11271 

z 

{^tr \ 

[2.150817 

[2.091612 

r)" j 

I 3.75927 
[3.660588 

A+ 1 

OL 1 

^1.871810 

1.580112 

El 

j 

f 32.24878 


f 34.05096 

J 

'30.37745 

1 

1 17.48596 


[13.38814 

1 

33.65919 

S2 

j 

1 32.20434 


1 34.00404 

J 

^30.33559 

1 

[17.57665 


[13.45758 

1 

^33.83377 

Xi 


f 11.02452 


[23.27915 


'33.31814 


[10.48005 


[2.853295 


31.94041 

El 


1-0.5247187 


[-0.31016527 


3.476590 


[-0.1283214 


[-2.547260 X 10“= 


0.2200185 

X2 


1 10.41513 


I 1.257832 


^9.400643 


[28.75634 


[11.51061 


7.017899 

E2 


[0.3775152 


[-8.204163 X 10“® 


-8.344217 X 


[-9.269701 X 10"^ 


[-0.2766008 


-9.616869 X 

^3 


[2.315848 


[8.313963 


33.06508 


[2.014284 


[30.25994 


0.2462918 

Ej 


1-1.307939 X 10-= 


[-0.1634056 


-3.258311 


[-6.897436 x 10"® 


[-0.1745266 


-7.002609 X 

X4 


I 39.95028 




— 


[53.07716 




J6.39366 

E4 


[-0.1030731 
[-3.027917 X 10"= 


— 


{1.508835 X 

Xe 





J 

'-4.048544 






-3.548035 
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TABLE II. Same as table ^ for the coexistence curve Ms, the susceptiblity x and the specific heat C calculated in the 
inhomogeneous phase (T < Tc and for n = 1). Notice that for each bound, the critical exponents 7 , a and the subcritical 
exponent A take on the same values as in table ^ (as it must according to the theory). An indication on the accuracy of the 
present work is provided by the parameter Xq (the critical background of the specific heat) which should take on the same 
value as in table | and differs slightly for the bound “max”. 



n = 1, inhomogeneous phase: A | 

“max”: 0.49862 
“min”: 0.50516 


Ms 



C 


e 


f 0.3270735 

1 0.3244954 

7 1 

f 1.2408875 
[1.23830 

—a < 

f -0.1049675 
[-0.11271 

z 


[0.93804691 

1^0.93700952 

{r-r \ 

1 18.386609 
[17.160196 

— < 

a 1 

1 3.3664988 
[3.0489086 

51 

j 

f 35.738988 

J 

f 2.3395295 


f 1.6036462 

i 

1 11.312578 

1 

[231.57911 


[4.7885026 

S2 

j 

1 35.689736 

J 

1 2.3363054 


1 1.6014362 

1 

[11.371253 

1 

[232.78025 


[4.8133395 


Xi 

r 303.21696 
[241.51662 

El 

1-1.7687565 x 10"' 
[-5.7056933 x lO"' 

X 2 

j 9.3779630 
[13.371447 

E 2 

j 0.17204103 
[0.20926322 

X 3 

J 1.3921229 
[248.39869 

Eg 

J6.0877711 X 10"® 
[-7.7226529 x 10"' 

X 4 

J 30.597947 
[82.917148 

E 4 

Jo. 17144092 
[0.15795660 

Xs 

J 6.5064180 
[4.7939978 

Es 

J-1.9479626 x 10 "' 
[4.8568967 x 10"^ 


1 76.549557 
23.010983 
-3.4865628 
0.88915951 
59.838911 
61.975912 
-16.395572 
4.0096724 
3.6904512 
316.29079 
4.7894058 x 10"^ 
-0.15361387 
63.029796 
50.582996 
19.215714 
-5.2595105 
9.3807398 
2.5909921 
0.13675156 
3.7692433 x 10"^ 


1 79.538017 
123.82899 
7.8643478 x 10"^ 
-0.26171236 
4.0631542 X 10"^ 
0.42099375 
9.1522424 x 10"® 
-3.6312569 x 10"® 
16.574905 
10.791900 
-0.28063252 
7.2072941 X 10"^ 
14.361662 
86.921949 
0.13070258 
0.41499782 
19.477188 
0.57982484 
0.28112994 
3.6928566 x 10"® 


f-4.0481532 
I-3.5480350 
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TABLE III. Same as table | for the correlation length the susceptiblity x a-nd the specific heat C calculated in the 
homogeneous phase (T > Tc and for n = 2). Comparing with tables ^ and one may notice the correlation of the value of 
the parameter Xe (the critical background of the specific heat) with the values of leading amplitude ja and of a: when a 
vanishes, and aX^ take on opposite values so as to transform the power law behavior |t|““ into the logarithmic singularity 
In |t|. 



n = 2, homogeneous phase: A j 

'“max”: 0.52551 
^ “min”: 0.52986 





C 


e 

V < 

[ 0.67181082 

1 0.66878932 

7 j 

[ 1.3188985 
[1.3148952 

—a < 

[ 1.5440 X 10"^ 
[6.370 X 10“® 



1 2.6289918 
[2.5496125 


1 5.5612909 
[5.3464216 

— \ 

a 1 

1-55.881907 

[-121.13056 

Si 


f 15.963748 

j 

f 96.831346 


f 4.0480920 


1 33.474847 

1 

[60.160224 


[28.884078 

S2 


1 16.381644 

j 

1 99.366178 


[4.1540622 


[34.505171 

1 

[62.011900 


[29.773102 

Xi 


f 28.529734 


f 16.310867 


f 1.1911602 


[111.52736 


[ 11.716728 


[37.837491 

Yi 


1-9.0963764 x 10“^ 


1-0.57358194 


1 5.5704053 x 10"^ 


[-2.8491431 X 10“^ 


1-4.1213479 X 10"^ 


[3.5434706 

X2 


1 9.1112497 


1 3.6615694 


1 1.2675164 


[13.427180 


[15.104245 


[59.951524 

Y2 


1-0.22311836 


1 5.6950360 x 10“^ 


[-5.8499557 x 10"' 


[-15.783043 


[-0.53630510 


[-1.7677012 X lO"' 

Xi 


1 0.11326011 


1 0.32669257 


1 27.562173 


[24.100833 


[3.1051883 


[1.3847300 

Yi 


1 3.8347877 x 10“^ 


1-4.8535318 x 10-“^ 


[-7.7243929 x 10"' 


[1.8612028 X 10“® 


[-3.8536260 x 10"^ 


[2.2769060 x 10"'‘ 

Xi 


1 72.907613 


[430.16727 


[46.806723 


[13.360107 


[239.95179 


[37.714984 

Yi 


1-5.0166740 X 10“^ 


1-6.7793463 x 10"® 


[-2.0360103 x 10~' 


[15.603262 


[-1.3735564 x 10"^ 


[-3.5387612 

Xi 


1 10.299474 
[7.3960558 


— 


— 

Y5 


1 2.0243740 x 10“^ 
[-0.13116817 


— 


— 

Xi 





J 

[50.158572 





1 

I 115.95104 
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TABLE IV. Same as table | for the correlation length the susceptiblity x s'-nd the specific heat C calculated in the 
homogeneous phase (T > Tc and for n = 3). 



n = 3, homogeneous phase: A j 

'“max”: 0.55227 
^ “min”: 0.55702 





C 


e 

V> < 

f 0.71090629 

1 0.70381062 

7 1 

f 1.3946000 
[1.3845100 

—a < 

f 0.132720 
[0.11143582 

z 

\ 

1 3.1722403 
[2.9632572 

in" \ 

1 7.9856105 
[7.2687650 

A+ 1 

a 1 

1-20.228436 

[-18.976690 

Si 

j 

f 31.477107 

J 

f 72.301387 


f 17.216858 

1 

1 78.590543 

1 

[52.048362 


[96640.818 

S2 

j 

1 33.213159 

J 

1 76.289014 


1 18.166416 

1 

[83.342746 

1 

[55.195616 


[ 102484.48 

Xi 


f 394.95293 


f 13.735280 


f 389.17897 


[10.931307 


[ 0.15890396 


[52.385639 

Yi 


1-5.2545625x10“® 


1-0.75390860 


1-3.0498190x10“ 


[16.025379 


[1.1525134x10“® 


[0.19945718 

X2 


1 0.15078920 


1 1.3616777 


1 4.7607464x10“^ 


[3.2123716 


[11.028433 


[1179.5468 

Y2 


1 2.8641031x10“® 


1 0.46960131 


1 3.2528813xl0“'‘ 


[-0.11355993 


[-4.7495792x10“^ 


[1.3062104x10“® 

X3 


1 11.387266 


1 1.3862969 


1 12.991689 


[505.59521 


[12.528570 


[22.441166 

Yi 


1-0.35982658 


1-0.49056204 


1-0.16565010 


[-6.9793726x10“® 


[-0.69601028 


[-0.13439992 

Xi 


1 78.089588 


1 437.65747 


1 65.185231 


[11.021201 


[312.81912 


[36.340773 

Yi 


1-5.6692560x10“^ 


1-1.4330672x10“^ 


1-0.10948246 


[-16.378555 


[-1.8816147x10“^ 


[-0.35127870 

Vs 


1 0.19582193 


1 - 


1 1.7080277 


[2.8188361 


[0.97748668 


[11.698778 

y5 


1-2.9029786x10“® 


1 - 


1 1.2417093x10“^ 


[6.6094313x10“^ 


[-7.8502947x10“® 


[6.2043595x10“^ 

Vs 





J 

f 8.2684338 





1 

[9.1558605 
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TABLE V. Above: Values of the (universal) critical exponents as they are accounted for by the crossover functions defined 
in tables The numbers given in parenthesis correspond to the GZ respective error-bound estimates. Below: the scaling 

relations structurally satisfied for each bound of the crossover functions defined in tables ^ - IV (the expected theoretical values 
are zero). Due to the practical necessity of using a small number of predefined criteria in the (unique) resummation method 
used, the scaling relations are (automatically) better satisfied in the present work than in the final upper and lower bounds 
of GZ (the numbers in parenthesis correspond to their bound estimates which have not particularly been determined so as to 
satisfy the scaling relations). 




critical exponent values 

n 

7 

V 

a. f} 

1 

1.2408875 (1.2409) 

0.631678 (0.6317) 

0.1049675 (0.105) 0.3270735 (0.3272) 

1.23830 (1.2383) 

0.6290975 (0.6291) 

0.11271 (0.113) 0.3244954 (0.3244) 


1.3188985 (1.3189) 

0.67181082 (0.6718) 

-0.01544 (-0.015) 


1.3148952 (1.3149) 

0.66878932 (0.6688) 

-0.00637 (-0.007) 


1.39460 (1.3945) 

0.71090629 (0.7108) 

-0.132720 (-0.132) 


1.38451 (1.3845) 

0.70381062 (0.7038) 

-0.11143582 (-0.112) 



scaling relations (should be zero) 


l«- 

-2 + 3i/| 

|a-f 2/3 -t 7 - 2 | 

1 

1.5 X 10’ 

(1. X 10"^) 

2. X 10"® (3. X 10"") 

2.5 X 10’ 

-® (3. X 10"^) 

8. X 10"^ (1. X IQ-*) 


7.5 X 10' 

■® (4. X 10"'^) 



2.0 X 10’ 

-® (6. X 10-^) 



1.1 X 10’ 

■® (4. X 10-^) 



4.0 X 10’ 

-® (6. X 10-"^) 



TABLE VI. Values of universal combinations of leading critical amplitudes for thermodynamic quantities combining calcu¬ 
lations in the two phases, hence for n = 1 only (from the crossover functions of tables and The two superposed numbers 
correspond to the respective bounds “max” (upper line) and “min” (lower line). In parenthesis are the respective bounds of 
the GZ estimates. 


A+/A- 

r+/r- 

R+ = A+V+jB'^ 

0.55601 (0.556) 

4.89100 (4.89) 

0.05940 (0.0594) 

0.51826 (0.518) 

4.68783 (4.69) 

0.05541 (0.0554) 


TABLE VII. Values of the universal quantity combi ning calculations in the homogeneous phase only, 

hence for three values of n (from the crossover functions of tables |, and p^ ) . The two superposed numbers correspond to 
the respective bounds “max” (upper line) and “min” (lower line). In parenthesis are rounded forms of the same estimates. 


n 

1 

2 

3 

Rt 

D 97n9Q 

0 : 2^99 (0-2696±0.0007) 

(0.3609±0.0012) 

O^SSe (0-4357±0.0025) 


18 





















TABLE VIII. Bounds on the fixed point values of the renormalized coupling g dehned as the zero of the Wilson function 
W{g). The resummation criteria for W{g) have been chosen so as to yield values of g* as close as possible to the GZ estimates 
given in parenthesis. The two superposed numbers correspond to the respective bounds “max” (upper line) and “min” (lower 
line). 


n 

1 

2 

3 

* 

1.41512 (1.415) 

1.40602 (1.406) 

1.39401 (1.394) 

9 

1.40687 (1.407) 

1.40004 (1.400) 

1.38605 (1.386) 


TABLE IX. Values of the universal ratios of the amplitudes of the first correction to scaling for thermodynamic quantities 
combining calculations in th e two phases, hence for n = 1 only (not obtained from the crossover functions, see text). Same 
presentation as in table |VIl[ 


/d-y 

^cl^c 

am 

(0.215±0.029) 

0.893 

1.823 (1-36±0.47) 

0 743 

0:048 


TABLE X. Values of the universal ratios of the amplitudes of the hrst correction to scaling for the quant ities calculated in 
the homogeneous phase (not obtained from the crossover functions, see text). Same presentation as in table 


VII. 


n 

1 

2 

3 

a+ /a+ <j 

(0 699 

0:659 

(0.636±0.006) 

a“ (»^«2 ±o,oi4) 

flj/aj 1 

(8.68±0.23) 

(5.97±0.13) 

(4.58±0.06) 
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APPENDIX A: COMPUTER PROGRAM IN FORTRAN 
FOR THE ISING MODEL IN THE TWO PHASES 
(UP-TO-DATE CRITERIA) 


program utdFTcrossl 

c up-to-date * 

c Crossover functions up to seven loop orders * 
c from Murray-Nickel (1991) and * 

c accounting for the recent analysis by Guida * 
c and Zinn-Justin [J Phys A31, 8103 (1998)] * 

c For Ising-like systems (n=l) above & below Tc * 

c see the companion paper: * 

C "Classical-to-critical crossovers from field * 
c theory" by C. Bagnuls and C. Bervillier * 

c The important subroutine is 
c fns(res,iphase,ibound,ifn,tau,theta,ikont) 
c res is the return value of the subroutine 
c iphase controls the phase: 
c iphase=l corresponds to T>Tc 

c iphase=2 corresponds to T<Tc 

c ifn controls the type of function chosen: 

c ifn=l gives the inverse correlation length for T>Tc if iphase=l 

c otherwise (iphase=2) it gives the coexistence curve 
c ifn=2 gives the inverse susceptibility 

c ifn=3 gives the specific heat 

c ibound controls the accuracy of the theoretical calculation: 
c ibound=l corresponds to "max" 

c ibound=2 corresponds to "min" 

c tau is IT-Tcl/Tc 

c theta is the adjustable parameter which relates the scaling field t to tau: 
c (in principle t=theta*tau but for practical use we have separated the asymptotic 

c pure scaling from the correction-to-scaling contributions, see the companion paper, 

c for theta=l one recovers the pure theoretical 

c functions of the scaling field t) 

c ikont is an integer which allows to choose the nature of res 
c if ikont=0 then res gives the functions as indicated above 

c if ikont=l then res corresponds to the effective exponent associated to the function 

c chosen according to the criteria described above, 

c All the functions are for Ising-like model (n=l) 

implicit double precision (a-h,o-z) 

C Evolution of the functions and of the effective exponents * 

C in terms of the pure scaling field t * 

C (i.e., theta=l) * 

C 

C Selection of the function 

iphase=2 ! controls the phase: l=T>Tc, 2=T<Tc 

ibound=l ! controls the accuracy bound : l=max, 2=min 

ifn=3 ! controls the fn: l=xi, 2=chi, 3=c; for T<Tc, l=coex curve 
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C End of selection 

C Files for saving the results 

open(29,F1LE=’expeff.dat’) 
open(30,FlLE=’Funct.dat’) 

theta=l ! gives the theoretical evolution in terms of the scaling field t 

t=l.d-9 

do 30 i=l,45 

call fns(f3,iphase,ibound,ifn,t,theta,0) !Call of the function with ikont=0 
write(30,666) t,f3 

call fns(exp2,iphase,ibound,ifn,t,theta,1) ! Call of the effective exponent with ikont=l 

print *, t,f3,exp2 
alogt=dloglO(t) 

write(29,666) alogt,exp2 ! Ready for a Log-Lin plot 

t=t*2 

30 continue 
close(29) 
close(30) 
stop 

666 format(2(G13.6,lx)) 

end 

C Theoretical functions for n=l 

SUBROUTINE fns(res,i,j,k,t,theta,ikont) 

c if ikont=0: res gives the function 
c if ikont=l: res gives the effective exponent 
c Table f contains the values of the functions’ parameters, 
c At the end of each set a comment allows to distinguish the 
c kind of function concerned. 

c For each set, the first value is the critical exponent value 
c (or minus this value for the specific heat), the second value 
c is the inverse of the critical amplitude for xi and chi, the 
c critical amplitude for aim and c, the third value is the critical 
c background which is non zero only for c. 

c The two values of the subcritical exponent delta corresponding to the two 
c bounds "max" and "min" are contained in the table del 

implicit double precision (a-h,o-z) 
dimension del(2),x(5),y(5) 
dimension f(15,3,2,2) 
data del/0.49862,0.50516/ 

data F/0.631678,2.150817,0.,32.24878,32.20434,11.02452, 


#-0.5247187,10.41513,0.3775152,2.315848,-1.3079390-02, 
#39.95028,-1.0307310-01,0.,0., 

!ximax 

T>Tc 

(corr length) 

#1.2408875,3.75927,0.,34.05096,34.00404,23.27915,-0.31016527, 

#1.257832,-8.2041630-03,8.313963,-0.1634056,0.,0.,0.,0., 

!chimax 

T>Tc 

(suscept) 

#-.1049675,1.871810,-4.048544,30.37745,30.33559,33.31814,3.476590, 
#9.400643,-8.3442170-03,33.06508,-3.258311,0.,0.,0.,0., 

!cmax T>Tc (specif heat) 

#0.6290975,2.091612,0.,17.48596,17.57665,10.48005,-0.1283214, 

#28.75634,-9.2697010-02,2.014284,-6.8974360-03,53.07716, 
#-3.0279170-02,0.,0., 

!ximin 

T>Tc 

(corr length) 

#1.23830,3.660588,0.,13.38814,13.45758,2.853295,-2.5472600-02, 
#11.51061,-0.2766008,30.25994,-0.1745266,0.,0.,0.,0., 

! chimin 

T>Tc 

(suscept) 
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#-.11271,1.580112,-3.548035,33.65919,33.83377,31.94041, 

#0.2200185,7.017899,-9.616869D-03,0.2462918,-7.002609D-05,76.39366, 
#1.5088350-02,0.,0., 

#0.32707350,0.93804691,0.,35.738988,35.689736, 

#303.21696,-0.176875650-02,9.3779630,0.17204103,1.3921229, 

#0.608777110-02,30.597947,0.17144092,6.5064180,-0.194796260-02, 

#1.2408875,18.386609,0.,2.3395295,2.3363054,76.549557,-3.4865628, 
#59.838911,-16.395572,3.6904512,0.478940580-01,63.029796, 

#19.215714,9.3807398,0.13675156, 

#-0.10496750,3.3664988,-4.0481532,1.6036462,1.6014362, 

#79.538017,0.786434780-01,0.406315420-01,0.915224240-04,16.574905, 
#-0.28063252,14.361662,0.13070258,19.477188,0.28112994, 

#0.32449540,0.93700952,0.,11.312578,11.371253, 

#241.51662,-0.570569330-01,13.371447,0.20926322,248.39869, 

#-0.772265290-02,82.917148,0.15795660,4.7939978,0.485689670-01, 

#1.2383000,17.160196,0.,231.57911,232.78025,23.010983, 
#0.88915951,61.975912,4.0096724,316.29079,-0.15361387, 

#50.582996,-5.2595105,2.5909921,0.376924330-01, 

#-0.11271000,3.0489086,-3.5480350,4.7885026,4.8133395, 

#123.82899,-0.26171236,0.42099375,-0.363125690-02,10.791900, 

#0.720729410-01,86.921949,0.41499782,0.57982484,0.369285660-02/ 

C definitions of constants 


d=del(j) ! correction exponent 

e=F(l,k,j,i) ! asymptotic exponent 

z=f(2,k,j,i) ! asymptotic amplitude 

x6=f(3,k,j,i) ! additive critical background 

sl=f(4,k,j,i) 

s2=f(5,k,j,i) 

kmax=0 

do 1 ii=l,5 

x(ii)=f(2*ii+4,k,j,i) ! definitions of the Xi’s and Yi’s 
if(x(ii).eq.0.) go to 2 
y(ii)=f(2*ii+5,k,j,i) 

1 kmax=kmax+l 

2 continue 

if(ikont.eq.0) then ! calculation of the function 
res=Z*t**e 


if(theta.eq.0.) then 

res=res+x6 

return 

end if 


If theta=0 

there is no correction-to-scaling 
only the pure scaling form 
survives 


tt=theta*t 

D=D-l+(Sl*dsqrt(tt)+l)/(S2*dsqrt(tt)+l) 

trr=(tt)**D 

do 3 kk=1,kmax 

3 res=res*(l+X(kk)*trr)**Y(kk) 
res=res+X6 


return 
end if 

if(ikont.eq.1) then 
if(theta.eq.0) then 
res=e 

if(k.eq.3) res=-res 
return 
end if 
tt=theta*t 


! calculation of the effective exponent 
If theta=0 
then 

the effective exponent 
reduces to the constant 
critical exponent 


! cmin T>Tc (specif heat) 

! aimax T<Tc (coex curve) 


! chimax T<Tc (suscept) 

! cmax T<Tc (specif heat) 

! aimin T<Tc (coex curve) 


! chimin T<Tc (suscept) 

! cmin T<Tc (specif heat) 
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D=D-l+(Sl*dsqrt(tt)+l)/(S2*dsqrt(tt)+l) 
trr=(tt)**D 

Dp=(S1-S2)/(2*dsqrt(tt)*(S2*dsqrt(tt)+l)**2) 
res=0. 

do 4 kk=1,kmax 

4 res=res+X(kk)*Y(kk)/(1+X(kk)*trr) 

res=res*(Dp*tt*DLOG(tt)+D)*trr 
res=res+e 

if(k.eq.3) res=-res ! the exponent changes sign for c 

return 

end if 

print *, "Error, ikont is ",ikont, ", but must be 0 or 1." 

stop 

end 
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APPENDIX B: COMPUTER PROGRAM IN FORTRAN 
FOR THE A- VECTOR MODEL (A = 1,2,3) IN THE HOMOGENEOUS PHASE 

(UP-TO-DATE CRITERIA) 


program utdFTcrossN 

c up-to-date * 

c Crossover functions up to seven loop orders * 

c from Murray-Nickel (1991) and * 

c accounting for the recent analysis by Guida * 

c and Zinn-Justin [J Phys A31, 8103 (1998)] * 

c For N-vector-like systems (n=l, 2, 3) * 

c in the homogeneous phase only * 

c see the companion paper: * 

C "Classical-to-critical crossovers from field * 

c theory" by C. Bagnuls and C. Bervillier * 

c The important subroutine is 
c fns(res,iN,ibound,ifn,tau,theta,ikont) 
c res is the return value of the subroutine 
c iN controls the value of n: 
c iN=l corresponds to n=l 

c iN=2 corresponds to n=2 

c iN=3 corresponds to n=3 

c ifn controls the type of function chosen: 
c ifn=l gives the inverse correlation length 

c ifn=2 gives the inverse susceptibility 

c ifn=3 gives the specific heat 

c ibound controls the accuracy of the theoretical calculation: 
c ibound=l corresponds to "max" 

c ibound=2 corresponds to "min" 

c tau is IT-Tcl/Tc 

c theta is the adjustable parameter which relates the scaling field t to tau: 
c (in principle t=theta*tau but for practical use we have separated the asymptotic 

c pure scaling from the correction-to-scaling contributions, see the companion paper, 

c for theta=l one recovers the pure theoretical 

c functions of the scaling field t) 

c ikont is an integer which allows to choose the nature of res 
c if ikont=0 then res gives the functions as indicated above 

c if ikont=l then res corresponds to the effective exponent associated to the function 

c chosen according to the criteria described above. 

implicit double precision (a-h,o-z) 

C Evolution of the functions and of the effective exponents * 

C in terms of the pure scaling field t * 

C (i.e., theta=l) * 

C 

C Selection of the function 

iN=2 ! controls the value of n: 1, 2, or 3 

ibound=l ! controls the accuracy bound : l=max, 2=min 

ifn=3 ! controls the fn: l=xi, 2=chi, 3=c 
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C End of selection 

C Files for saving the results 

open(29,F1LE=’expeff.dat’) 
open(30,FlLE=’Funct.dat’) 

theta=l ! gives the theoretical evolution in terms of the scaling field t 

t=l.d-9 

do 30 i=l,45 

call fns(f3,iN,ibound,ifn,t,theta,0) !Call of the function with ikont=0 
write(30,666) t,f3 

call fns(exp2,iN,ibound,ifn,t,theta,1) ! Call of the effective exponent with ikont=l 

print *, t,f3,exp2 
alogt=dloglO(t) 

write(29,666) alogt,exp2 ! Ready for a Log-Lin plot 

t=t*2 

30 continue 
close(29) 
close(30) 
stop 

666 format(2(G13.6,lx)) 

end 

C Theoretical functions for n-vector-like systems 

SUBROUTINE fns(res,i,j,k,t,theta,ikont) 

c if ikont=0: res gives the function 
c if ikont=l: res gives the effective exponent 
c Table f contains the values of the functions’ parameters, 
c At the end of each set a comment allows to distinguish the 
c kind of function concerned. 

c For each set, the first value is the critical exponent value 
c (or minus this value for the specific heat), the second value 
c is the inverse of the critical amplitude for xi and chi, the 
c critical amplitude for c, the third value is the critical 
c background which is non zero only for c. 

c The six values of the subcritical exponent delta corresponding to the two 
c bounds "max" and "min" for each value of n are contained in the table del 

implicit double precision (a-h,o-z) 
dimension del(2,3),x(5),y(5) 
dimension f(15,3,2,3) 

data del/0.49862,0.50516,0.52551,0.52986,0.55227,0.55702/ 
data F/0.631678,2.150817,0.,32.24878,32.20434,11.02452, 
#-0.5247187,10.41513,0.3775152,2.315848,-1.3079390-02, 

#39.95028,-1.0307310-01,0.,0., 

#1.2408875,3.75927,0.,34.05096,34.00404,23.27915,-0.31016527, 

#1.257832,-8.2041630-03,8.313963,-0.1634056,0.,0.,0.,0., 

#-.1049675,1.871810,-4.048544,30.37745,30.33559,33.31814,3.476590, 
#9.400643,-8.3442170-03,33.06508,-3.258311,0.,0.,0.,0., 

#0.6290975,2.091612,0.,17.48596,17.57665,10.48005,-0.1283214, 

#28.75634,-9.2697010-02,2.014284,-6.8974360-03,53.07716, 

#-3.0279170-02,0.,0., 

#1.23830,3.660588,0.,13.38814,13.45758,2.853295,-2.5472600-02, 
#11.51061,-0.2766008,30.25994,-0.1745266,0.,0.,0.,0., 


Iximax n=l (corr length) 

! chimax n=l (suscept) 

! cmax n=l (specif heat) 

Iximin n=l (corr length) 
! chimin n=l (suscept) 
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#-.11271,1.580112,-3.548035,33.65919,33.83377,31.94041, 

#0.2200185,7.017899,-9.616869D-03,0.2462918,-7.002609D-05,76.39366, 
#1.5088350-02,0.,0., 

#0.67181082,2.6289918,0.,15.963748,16.381644, 

#28.529734,-0.909637640-01,9.1112497,-0.22311836,0.11326011, 

#0.383478770-03,72.907613,-0.501667400-01,10.299474,0.202437400-01, 
#1.3188985,5.5612909,0.,96.831346,99.366178, 
#16.310867,-0.57358194,3.6615694,-0.569503600-01,0.32669257, 

#-0.485353180-03,430.16727,-0.677934630-02,0 . , 0. , 

#0.154400000-01,-55.881907,50.158572,4.0480920,4.1540622, 

#1.1911602,0.557040530-01,1.2675164,-0.584995570-01,27.562173, 

#-0.772439290-02,46.806723,-0.203601030-01,0.,0., 

#0.66878932,2.5496125,0.,33.474847,34.505171, 

#111.52736,-0.284914310-01,13.427180,-15.783043,24.100833, 

#0.186120280-02,13.360107,15.603262,7.3960558,-0.13116817, 

#1.3148952,5.3464216,0.,60.160224,62.011900, 

#11.716728,-0.412134790-01,15.104245,-0.53630510,3.1051883, 

#-0.385362600-01,239.95179,-0.137355640-01,0.,0. , 

#0.637000000-02,-121.13056,115.95104,28.884078,29.773102, 

#37.837491,3.5434706,59.951524,-0.176770120-01,1.3847300, 
#0.227690600-03,37.714984,-3.5387612,0.,0., 

#0.71090629,3.1722403,0.,31.477107,33.213159, 

#394.95293,-0.525456250-02,0.15078920,0.286410310-02,11.387266, 

#-0.35982658,78.089588,-0.566925600-01,0.19582193,-0.290297860-02, 
#1.3946000,7.9856105,0.,72.301387,76.289014, 

#13.735280,-0.75390860,1.3616777,0.46960131,1.3862969, 
#-0.49056204,437.65747,-0.143306720-01,0.,0., 

#0.13272000,-20.228436,8.2684338,17.216858,18.166416, 

#389.17897,-0.304981900-02,0.476074640-01,0.325288130-03,12.991689, 
#-0.16565010,65.185231,-0.10948246,1.7080277,0.124170930-01, 

#0.70381062,2.9632572,0.,78.590543,83.342746, 

#10.931307,16.025379,3.2123716,-0.11355993,505.59521, 

#-0.697937260-02,11.021201,-16.378555,2.8188361,0.660943130-01, 

#1.3845100,7.2687650,0.,52.048362,55.195616, 

#0.15890396,0.115251340-02,11.028433,-0.474957920-01,12.528570, 

#-0.69601028,312.81912,-0.188161470-01,0.97748668,-0.785029470-02, 
#0.11143582,-18.976690,9.1558605,96640.818,102484.48, 

#52.385639,0.19945718,1179.5468,0.130621040-02,22.441166, 

#-0.13439992,36.340773,-0.35127870,11.698778,0.620435950-01/ 

C definitions of constants 


if(i.gt.3.or.i.It.1) then 

print *, "Error, n= ",i, "is not available. Sorry. (0<n<4)" 

stop 

end if 

d=del(j,i) ! correction exponent 

e=F(l,k,j,i) ! asymptotic exponent 

z=f(2,k,j,i) ! asymptotic amplitude 

x6=f(3,k,j,i) ! additive critical background 

sl=f(4,k,j,i) 

s2=f(5,k,j,i) 

kmax=0 

do 1 ii=l,5 

x(ii)=f(2*ii+4,k,j,i) ! definitions of the Xi’s and Yi’s 
if(x(ii).eq.0.) go to 2 
y(ii)=f(2*ii+5,k,j,i) 

1 kmax=kmax+l 


! cmin n- 


! ximax n- 


!chimax 


! cmax n=: 


! ximin 


! chimin 


! cmin n- 


! ximax n- 


!chimax 


! cmax n=: 


! ximin 


! chimin 


! cmin n- 


■1 (specif heat) 
•2 (corr length) 
=2 (suscept) 

I (specif heat) 


1=2 (corr length) 
n=2 (suscept) 

•2 (specif heat) 


•3 (corr length) 
l= 3 (suscept) 

1 (specif heat) 

=3 (corr length) 


n=3 (suscept) 
•3 (suscept) 
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2 continue 

if(ikont.eq.0) then ! calculation of the function 
res=Z*t**e 

if(theta.eq.0.) then ! If theta=0 

res=res+x6 ! there is no correction-to-scaling 

return ! only the pure scaling form 

end if ! survives 

tt=theta*t 

D=D-l+(Sl*dsqrt(tt)+l)/(S2*dsqrt(tt)+l) 

trr=(tt)**D 

do 3 kk=1,kmax 

3 res=res*(l+X(kk)*trr)**Y(kk) 
res=res+X6 

return 
end if 

if(ikont.eq.1) then ! calculation of the effective exponent 
if(theta.eq.0) then ! If theta=0 

res=e ! then 

if(k.eq.3) res=-res ! the effective exponent 
return ! reduces to the constant 

end if ! critical exponent 

tt=theta*t 

D=D-l+(Sl*dsqrt(tt)+l)/(S2*dsqrt(tt)+l) 
trr=(tt)**D 

Dp=(S1-S2)/(2*dsqrt(tt)*(S2*dsqrt(tt)+l)**2) 
res=0. 

do 4 kk=1,kmax 

4 res=res+X(kk)*Y(kk)/(1+X(kk)*trr) 
res=res*(Dp*tt*DL0G(tt)+D)*trr 
res=res+e 

if(k.eq.3) res=-res ! the exponent changes sign for c 

return 

end if 

print *, "Error, ikont is ",ikont, ", but must be 0 or 1." 

stop 

end 
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APPENDIX C: COMPUTER PROGRAM IN FORTRAN 
FOR THE ISING MODEL IN THE TWO PHASES 
(FORMER CRITERIA, CORRECTED SERIES) 


program FTcrossl 


c former corrected * 
c Crossover functions up to six loop orders * 
c after correction of the errors for the * 
c inhomogeneous phase * 


c For Ising-like systems (n=l) above & below Tc * 

c see the companion paper: * 

C "Classical-to-critical crossovers from field * 
c theory" by C. Bagnuls and C. Bervillier * 

c The important subroutine is 
c fns(res,iphase,ibound,ifn,tau,theta,ikont) 
c res is the return value of the subroutine 
c iphase controls the phase: 
c iphase=l corresponds to T>Tc 

c iphase=2 corresponds to T<Tc 

c ifn controls the type of function chosen: 

c ifn=l gives the inverse correlation length for T>Tc if iphase=l 

c otherwise (iphase=2) it gives the coexistence curve 
c ifn=2 gives the inverse susceptibility 

c ifn=3 gives the specific heat 

c ibound controls the accuracy of the theoretical calculation: 
c ibound=l corresponds to "max" 

c ibound=2 corresponds to "min" 

c tau is IT-Tcl/Tc 

c theta is the adjustable parameter which relates the scaling field t to tau: 
c (in principle t=theta*tau but for practical use we have separated the asymptotic 

c pure scaling from the correction-to-scaling contributions, see the companion paper, 

c for theta=l one recovers the pure theoretical 

c functions of the scaling field t) 

c ikont is an integer which allows to choose the nature of res 
c if ikont=0 then res gives the functions as indicated above 

c if ikont=l then res corresponds to the effective exponent associated to the function 

c chosen according to the criteria described above, 

c All the functions are for Ising-like model (n=l) 

implicit double precision (a-h,o-z) 

C Evolution of the functions and of the effective exponents * 

C in terms of the pure scaling field t * 

C (i.e., theta=l) * 

C 

C Selection of the function 

iphase=2 ! controls the phase: l=T>Tc, 2=T<Tc 

ibound=l ! controls the accuracy bound : l=max, 2=min 

ifn=3 ! controls the fn: l=xi, 2=chi, 3=c; for T<Tc, l=coex curve 
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C End of selection 

C Files for saving the results 

open(29,F1LE=’expeff.dat’) 
open(30,FlLE=’Funct.dat’) 

theta=l ! gives the theoretical evolution in terms of the scaling field t 

t=l.d-9 

do 30 i=l,45 

call fns(f3,iphase,ibound,ifn,t,theta,0) !Call of the function with ikont=0 
write(30,666) t,f3 

call fns(exp2,iphase,ibound,ifn,t,theta,1) ! Call of the effective exponent with ikont=l 

print *, t,f3,exp2 
alogt=dloglO(t) 

write(29,666) alogt,exp2 ! Ready for a Log-Lin plot 

t=t*2 

30 continue 
close(29) 
close(30) 
stop 

666 format(2(G13.6,lx)) 

end 

C Theoretical functions for n=l 

SUBROUTINE fns(res,i,j,k,t,theta,ikont) 

c if ikont=0: res gives the function 
c if ikont=l: res gives the effective exponent 
c Table f contains the values of the functions’ parameters, 
c At the end of each set a comment allows to distinguish the 
c kind of function concerned. 

c For each set, the first value is the critical exponent value 
c (or minus this value for the specific heat), the second value 
c is the inverse of the critical amplitude for xi and chi, the 
c critical amplitude for aim and c, the third value is the critical 
c background which is non zero only for c. 

c The two values of the subcritical exponent delta corresponding to the two 
c bounds "max" and "min" are contained in the table del 

implicit double precision (a-h,o-z) 
dimension del(2),x(5),y(5) 
dimension f(15,3,2,2) 
data del/0.49125,0.50031/ 

data F/0.630501,2.12411,0.,16.87409,16.72772,23.51302,-0.2154827, 

#0.6009808,-1.955281D-03,5.460931,-4.356399D-02,0.,0.,0.,0., 

#1.24194,3.80403,0.,10.7784,10.6849,18.2831,-0.447952,0.211862, 

#-l.12883D-03,2.82005,-3.47989D-02,0.,0.,0.,0., 

#-.108496,1.74928,-3.87997,22077.0,21885.5,8.15287,-8.12340D-03, 
#36.4432,0.223350,395.406,1.765200-03,0.,0.,0.,0., 

#0.629121,2.09256,0.,1.793553,1.794109,32.37368,-0.1214804, 
#11.42729,-0.1263489,2.192044,-1.0412700-02,0.,0.,0. ,0. , 

#1.239485,3.706359,0.,1.229767,1.230149,2.198050,-2.0154540-02, 
#24.55992,-0.2546755,10.44460,-0.2041400,0.,0.,0.,0., 

#-.112636,1.58819,-3.57607,69.75599,69.77762,11.38361, 

#-l.0049830-01,11.37558,8.0309140-02,33.34772,0.2454611,0.,0.,0., 

#0. , 


Iximax T>Tc (corr length) 
! chimax T>Tc (suscept) 

! cmax T>Tc (specif heat) 
Iximin T>Tc (corr length) 
! chimin T>Tc (suscept) 

! cmin T>Tc (specif heat) 
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#0.32516,0.9213251,0.,64.93122,64.36800,1.770551,1.028455D-02, 

#46.28114,0.1928966,89.74606,-5.501773D-02,10.80746,0.2015166,0., 
#0. , 

#1.24194,19.19139,0.,8.662018,8.586883,8.293155,0.2623473, 

#37.19002,-1.070686,71.02634,-2.254137,62.54640,2.578595, 

#0. ,0. , 

#-.108496,3.184180,-3.879732,49.56056,49.13067,4.262502, 

#6.1456000-03,66.44198,0.1076623,20.11234,1.0318410-01,0.,0.,0.,0., 
#0.32357,0.8971946,0.,13.29911,13.30323,131.5962,2.4309340-02, 

#1.304663,4.7346440-03,23.57031,0.1913076,8.270542,0.1325084,0. , 0 . , 
#1.239485,16.78313,0.,360.4357,360.5474,22.85732,0.4525390, 

#204.6138,5.7880260-02,3.581834,6.5275020-02,41.70674,-1.054664, 

#0. ,0. , 

#-0.112636,3.004344,-3.576325,225.67103,225.74101,13.429043, 

#7.6242080-02,183.93336,-29.012144,56.400856,0.16278169,183.93769, 
#28.998392,0.,0./ 

C definitions of constants 

d=del(j) ! correction exponent 

e=F(l,k,j,i) ! asymptotic exponent 

z=f(2,k,j,i) ! asymptotic amplitude 

x6=f(3,k,j,i) ! additive critical background 

sl=f(4,k,j,i) 

s2=f(5,k,j,i) 

kmax=0 

do 1 ii=l,5 

x(ii)=f(2*ii+4,k,j,i) ! definitions of the Xi’s and Yi’s 
if(x(ii).eq.0.) go to 2 
y(ii)=f(2*ii+5,k,j,i) 

1 kmax=kmax+l 

2 continue 

if(ikont.eq.0) then ! calculation of the function 
res=Z*t**e 

if(theta.eq.0.) then ! If theta=0 

res=res+x6 ! there is no correction-to-scaling 

return ! only the pure scaling form 

end if ! survives 

tt=theta*t 

D=D-l+(Sl*dsqrt(tt)+l)/(S2*dsqrt(tt)+l) 

trr=(tt)**D 

do 3 kk=1,kmax 

3 res=res*(l+X(kk)*trr)**Y(kk) 
res=res+X6 

return 
end if 

if(ikont.eq.1) then ! calculation of the effective exponent 
if(theta.eq.0) then ! If theta=0 

res=e ! then 

if(k.eq.3) res=-res ! the effective exponent 

return ! reduces to the constant 

end if ! critical exponent 

tt=theta*t 

D=D-l+(Sl*dsqrt(tt)+l)/(S2*dsqrt(tt)+l) 
trr=(tt)**D 

Dp=(Sl-S2)/(2*dsqrt(tt)*(S2*dsqrt(tt)+l)**2) 
res=0. 

do 4 kk=1,kmax 


! aimax T<Tc (coex curve) 

! chimax T<Tc (suscept) 

! cmax T<Tc (specif heat) 
! aimin T<Tc (coex curve) 

! chimin T<Tc (suscept) 

! cmin T<Tc (specif heat) 
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4 


res=res+X(kk)*Y(kk)/(1+X(kk)*trr) 
res=res*(Dp*tt*DLOG(tt)+D)*trr 
res=res+e 

if(k.eq.3) res=-res ! the exponent changes sign for c 

return 

end if 

print *, "Error, ikont is ",ikont, ", but must be 0 or 1." 

stop 

end 
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APPENDIX D: COMPUTER PROGRAM IN FORTRAN 
FOR THE A- VECTOR MODEL (A = 1,2,3) IN THE HOMOGENEOUS PHASE 

(FORMER CRITERIA) 


progreun FTcrossN 

c former * 

c Crossover functions up to six loop orders * 

c For N-vector-like systems (n=l, 2, 3) * 

c in the homogeneous phase only * 

c see the companion paper: * 

C "Classical-to-critical crossovers from field * 
c theory" by C. Bagnuls and C. Bervillier * 

c The important subroutine is 
c fns(res,iN,ibound,ifn,tau,theta,ikont) 
c res is the return value of the subroutine 
c iN controls the value of n: 
c iN=l corresponds to n=l 

c iN=2 corresponds to n=2 

c iN=3 corresponds to n=3 

c ifn controls the type of function chosen: 
c ifn=l gives the inverse correlation length 

c ifn=2 gives the inverse susceptibility 

c ifn=3 gives the specific heat 

c ibound controls the accuracy of the theoretical calculation: 
c ibound=l corresponds to "max" 

c ibound=2 corresponds to "min" 

c tau is IT-Tcl/Tc 

c theta is the adjustable parameter which relates the scaling field t to tau: 
c (in principle t=theta*tau but for practical use we have separated the asymptotic 

c pure scaling from the correction-to-scaling contributions, see the companion paper, 

c for theta=l one recovers the pure theoretical 

c functions of the scaling field t) 

c ikont is an integer which allows to choose the nature of res 
c if ikont=0 then res gives the functions as indicated above 

c if ikont=l then res corresponds to the effective exponent associated to the function 

c chosen according to the criteria described above. 

implicit double precision (a-h,o-z) 

C Evolution of the functions and of the effective exponents * 

C in terms of the pure scaling field t * 

C (i.e. , theta=l) * 

C 

C Selection of the function 


iN=2 

! controls 

the 

value of 

n: 1, 2 

or 3 

ibound=l 

! controls 

the 

accuracy 

bound : 

l=max, 2=min 

ifn=3 

! controls 

the 

fn: l=xi 

2=chi, 

3=c 


C End of selection 


32 




C Files for saving the results 

open(29,F1LE=’expeff.dat’) 
open(30,FlLE=’Funct.dat’) 

theta=l ! gives the theoretical evolution in terms of the scaling field t 

t=l.d-9 

do 30 i=l,45 

call fns(f3,iN,ibound,ifn,t,theta,0) !Call of the function with ikont=0 
write(30,666) t,f3 

call fns(exp2,iN,ibound,ifn,t,theta,1) ! Call of the effective exponent with ikont=l 

print *, t,f3,exp2 
alogt=dloglO(t) 

write(29,666) alogt,exp2 ! Ready for a Log-Lin plot 

t=t*2 

30 continue 
close(29) 
close(30) 
stop 

666 format(2(G13.6,lx)) 

end 

C Theoretical functions for n-vector-like systems 

SUBROUTINE fns(res,i,j,k,t,theta,ikont) 

c if ikont=0: res gives the function 
c if ikont=l: res gives the effective exponent 
c Table f contains the values of the functions’ parameters, 
c At the end of each set a comment allows to distinguish the 
c kind of function concerned. 

c For each set, the first value is the critical exponent value 
c (or minus this value for the specific heat), the second value 
c is the inverse of the critical amplitude for xi and chi, the 
c critical amplitude for c, the third value is the critical 
c background which is non zero only for c. 

c The six values of the subcritical exponent delta corresponding to the two 
c bounds "max" and "min" for each value of n are contained in the table del 

implicit double precision (a-h,o-z) 
dimension del(2,3),x(5),y(5) 
dimension f(15,3,2,3) 

data del/0.49125,0.50031,0.52012,0.52737,0.54975,0.55043/ 
data F/0.630501,2.12411,0.,16.87409,16.72772,23.51302,-0.2154827, 

#0.6009808,-1.955281D-03,5.460931,-4.356399D-02,0.,0.,0.,0., 

#1.24194,3.80403,0.,10.7784,10.6849,18.2831,-0.447952,0.211862, 

#-l.12883D-03,2.82005,-3.47989D-02,0.,0.,0.,0., 

#-.108496,1.74928,-3.87997,22077.0,21885.5,8.15287,-8.12340D-03, 
#36.4432,0.223350,395.406,1.765200-03,0.,0.,0.,0., 

#0.629121,2.09256,0.,1.793553,1.794109,32.37368,-0.1214804, 
#11.42729,-0.1263489,2.192044,-1.0412700-02,0.,0.,0. ,0. , 

#1.239485,3.706359,0.,1.229767,1.230149,2.198050,-2.0154540-02, 
#24.55992,-0.2546755,10.44460,-0.2041400,0.,0.,0.,0., 

#-.112636,1.58819,-3.57607,69.75599,69.77762,11.38361, 

#-l.0049830-01,11.37558,8.0309140-02,33.34772,0.2454611,0.,0.,0., 

#0. , 

#0.669848,2.578296,0.,28.84310,29.43534,15.75544,-0.2585115, 

#5.012764,-4.3432500-02,77.73647,-3.7751960-02,0.,0.,0.,0., 


Iximax T>Tc n=l (corr length) 

! chimax T>Tc n=l (suscept) 

! cmax T>Tc n=l (specif heat) 
Iximin T>Tc n=l (corr length) 
! chimin T>Tc n=l (suscept) 

! cmin T>Tc n=l (specif heat) 
Iximax T>Tc n=2 (corr length) 


33 



#1.31792,5.51529,0.,63.8357,65.1464,16.9739,-0.508216, 

#1.85577,-1.44706D-02,7.33240,-1.03476D-01,247.042,-9.67699D-03, 
#0. ,0. , 

#9.488051D-03,-85.85105,80.36121,457.8256,467.2262,21.93721, 

#-4.630859D-02,11.62941,-9.686347D-03,15.18151,3.701884D-02, 

#0.,0.,0.,0. , 

#0.6678685,2.527959,0.,92.02055,94.61002,1.400964,-3.366864D-03, 
#5.520135,-4.367003D-02,395.4183,-4.673720D-03,18.17701, 
#-0.2840264,0.,0., 

#1.314001,5.307447,0.,60.60521,62.31065,15.15227,-0.5613569, 

#2.473990,-2.203622D-02,5.819209,-3.266531D-02,246.3158, 
#-1.1943530-02,0.,0., 

#3.6032330-03,-207.0019,201.9794,367.4423,377.7822,17.09462, 
#-12.10980,3.267621,-1.6635250-04,17.08843,12.10276,0.,0.,0.,0., 
#0.7056167,3.014192,0.,76.15010,80.13691,4.745939,-8.7775040-02, 
#445.7242,-7.1869120-03,8.077017,0.2125871,13.25811,-0.5288586, 
#0. ,0. , 

#1.389161,7.591932,0.,27.84825,29.30624,10.79849,-0.6921719, 

#5.6066450-02,5.5195550-04,316.7178,-8.4499680-03,67.73494, 
#-7.8251080-02,0.,0., 

#0.1168480,-19.46045,8.971331,69.91117,73.57134,389.0545, 

#-7.7599790-03,15.08428,-1.794460,0.2445648,3.6241320-04, 
#14.03141,1.568162,0.,0., 

#0.7038317,2.969913,0.,78.53203,82.70273,12.80904,-0.5435931, 
#5.710801,-0.3817143,463.4073,-6.8429270-03,6.792342,0.5244869, 
#0. ,0. , 

#1.383616,7.231499,0.,48.19012,50.74942,12.44953,-0.7397974, 

#0.5892605,0.2315015,0.6046381,-0.2390366,260.9201, 
#-1.9900060-02,0.,0., 

#0.1114902,-19.09212,9.185812,2.689863,2.832717,40.48173, 

#-0.1850101,1.173960,-0.1448486,1.7434640-02,6.1136890-04, 
#0.8190587,0.1062669,0.,0./ 


C definitions of constants 

if(i.gt.3.or.i.It.1) then 

print *, "Error, n= ",i, "is not available. Sorry. (0<n<4)" 

stop 

end if 


d=del(j,i) ! correction exponent 

e=F(l,k,j,i) ! asymptotic exponent 

z=f(2,k,j,i) ! asymptotic amplitude 

x6=f(3,k,j,i) ! additive critical background 

sl=f(4,k,j,i) 

s2=f(5,k,j,i) 

kmax=0 


1 

2 


do 1 ii=l,5 

x(ii)=f(2*ii+4,k,j ,i) 

if(x(ii).eq.0.) go to 2 

y(ii)=f(2*ii+5,k,j,i) 

kmax=kmax+l 

continue 

if(ikont.eq.0) then 

res=Z*t**e 

if(theta.eq.0.) then 

res=res+x6 

return 

end if 


! definitions of the Xi’s and Yi’s 


calculation of the function 
If theta=0 

there is no correction-to-scaling 
only the pure scaling form 
survives 


! chimax T>Tc n=2 (suscept) 

! cmax T>Tc n=2 (specif heat) 

Iximin T>Tc n=2 (corr length) 

! chimin T>Tc n=2 (suscept) 

! cmin T>Tc n=2 (specif heat) 

Iximax T>Tc n=3 (corr length) 

! chimax T>Tc n=3 (suscept) 

! cmax T>Tc n=3 (specif heat) 

Iximin T>Tc n=3 (corr length) 

! chimin T>Tc n=3 (suscept) 

! cmin T>Tc n=3 (suscept) 
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tt=theta*t 

D=D-l+(Sl*dsqrt(tt)+l)/(S2*dsqrt(tt)+l) 

trr=(tt)**D 

do 3 kk=1,kmax 

3 res=res*(l+X(kk)*trr)**Y(kk) 
res=res+X6 

return 
end if 

if(ikont.eq.1) then ! calculation of the effective exponent 
if(theta.eq.0) then ! If theta=0 

res=e ! then 

if(k.eq.3) res=-res ! the effective exponent 

return ! reduces to the constant 

end if ! critical exponent 

tt=theta*t 

D=D-l+(Sl*dsqrt(tt)+l)/(S2*dsqrt(tt)+l) 
trr=(tt)**D 

Dp=(S1-S2)/(2*dsqrt(tt)*(S2*dsqrt(tt)+l)**2) 
res=0. 

do 4 kk=1,kmax 

4 res=res+X(kk)*Y(kk)/(1+X(kk)*trr) 
res=res*(Dp*tt*DL0G(tt)+D)*trr 
res=res+e 

if(k.eq.3) res=-res ! the exponent changes sign for c 

return 

end if 

print *, "Error, ikont is ",ikont, ", but must be 0 or 1." 

stop 

end 
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